Theory of isotope effect in photoemission spectra of high-T c superconducting cuprates 
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We investigate the effect of isotope substitution on the eiectronic spectral functions within a 
model where the charge carriers are coupled to bosonic charge-order (CO) fluctuations centered 
around some mean frequency wo and with enhanced scattering at wave- vector q c . It is shown that 
a mass dependence of uio is not sufficient in order to account, especially at high energies, for the 
dispersion shifts experimentally observed in an optimally doped superconducting cuprate. We argue 
that isotope substitution induces a change of the spatial CO correlations which gives good account 
of the experimental data. 

PACS numbers: 74.72.-h, 74.20.Mn, 71.38.Cn, 79.60.-i 



Several isotope effects (IE) have been found in under- 
doped superconducting cuprates, not only in the critical 
temperature but also in the superfluid density, in the 
fermion effective mass, and in the dynamic lattice distor- 
tion and in the pseudogap crossover temperature T* 
Q ■ More recently a remarkable IE has also been found in 
angle-resolved photoemission spectroscopy (ARPES) ex- 
periments on optimally doped Bi-2212 |3j. In this exper- 
iment the dispersive single-particle excitations were in- 
vestigated and strong peculiar variations upon replacing 
16 O with 18 O were found. Specifically, little isotopic de- 
pendence was found both for the kink and the low-energy 
part of the spectra along the r — > X, [i.e., (0, 0) — > (w, tt)] 
direction. On the other hand the dispersive "hump" 
at higher binding energy displayed a sizable IE, which 
caused a shift of about 15 meV towards larger binding 
energies. Around the M points at (tt, 0) and (0, tt) the IE 
has a reversed sign and both the low- and high-energy 
parts of the dispersion shift to smaller binding energies 
upon isotope substitution. In this region shifts as large 
as 40 meV are reported. This IE is absolutely non triv- 
ial in various respects. First of all, the ion replacement 
is expected to change by only 6% the phonon frequen- 
cies involving oxygen ions. The changes of the electron 
dispersions are large and within a standard perturbative 
Eliashberg-like scheme would imply too large e-ph cou- 
plings 0| . Secondly the isotopic shifts have different signs 
depending on the region of the Fermi surface. Thirdly 
they do not substantially occur around the "kinks" of 
the dispersion curves, as it would be expected if the iso- 
tope replacement were only changing the frequency of the 
phonon [or more generally of the collective mode (CM)] 
interacting with the quasiparticles (QP's). 

In this letter we specifically address the issue of this 
peculiar IE in ARPES and we provide a natural expla- 
nation for it within the theory based on the occurrence 
of a charge-ordering (CO) quantum-critical point (QCP) 
around optimal doping of the cuprates 0. This theory 
predicts the occurrence of a CO phase in a substantial 



region of the phase diagram below a doping-dependent 
critical temperature Tco{x) ending at zero tempera- 
ture in a QCP around optimal doping. Actually sev- 
eral mechanisms (pairing, disorder, and low dimension- 
ality) may prevent a fully realized long-range CO state. 
Nevertheless, despite this "missed" criticality, in the low- 
temperature underdoped region below T* ~ TqOi local 
dynamical CO should be revealed as a natural inclination 
of the system. After this theoretical prediction, many ev- 
idences have arisen that a QCP could be hidden under 
the superconducting "dome" 0, 0] . The prediction of a 
local CO in the underdoped cuprates is also finding new 
experimental confirmations Q , which corroborate less re- 
cent observations of "stripe" phases. Within the CO- 
QCP theory it was also possible to explain the peculiar 
IE of the pseudogap crossover temperature T* mentioned 
above [[J . Therefore it seems legitimate to adopt this the- 
oretical framework to address the issue of IE in ARPES 
experiments. In particular there are several characteris- 
tics of the theory, which render it a particularly suitable 
scheme: a) the presence of a "missed" criticality justi- 
fies the presence of low-energy (nearly critical) charge 
collective modes (CM), which have already been shown 
to produce "kinks" in the dispersive spectra |l0| : b) on 
general grounds, the proximity to a ("missed") criticality 
renders the system highly susceptible so that large effects 
(e.g. large IE on the single-particle excitations) can arise 
from small changes (like a 6% shift of some phonon fre- 
quencies); c) within the Hubbard- Holstcin model, which 
was the starting point of the CO-QCP theory, the CM are 
mixed phonon-plasmon modes, which therefore have both 
a built-in dependence on the phonon frequency naturally 
leading to IE's and an electronic character, which renders 
them similar to the spin modes considered in Refs.[ll|: 
d) the CM's are only critical near specific wavevectors 
q c related to the wavelengths of the CO textures. As a 
consequence they have a strong momentum dependence, 
which allows them to play different roles in interacting 
with QP's in different regions of the Brillouin zone. As 
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we shall see this is a crucial feature in order to explain 
the opposite sign of the IE on dispersive bands along 
the nodal and antinodal directions. Within the above 
framework we performed a systematic analysis of the ef- 
fects of changing the CM frequency u)q, the QP-CM cou- 
pling g, and the CO inverse correlation length 7. Our 
main finding is that changes of luq and/or g alone by 
no means allow to reproduce the experimental results. 
On the other hand the crucial parameter to vary in or- 
der to semiquantitatively reproduce the experimental iso- 
topic shift is the inverse charge-charge coherence length 
7. Before reporting on our detailed calculations, we first 
present qualitative arguments to give an intuitive under- 
standing of this result. We start, e.g, from an hole state 
at wavevector p, which is scattered by a narrow CM of 
energy ujq and wavevector q c . This hole acquires a self- 
energy Ep(cj) = g 2 /(u + ujq — £ P +q c ). If the CM has 
instead a wavevector broadening 7, the self-energy be- 
comes proportional to the hole density of states (DOS) 
p p+qc in the region of the momenta, where the hole is 
scattered by the CM. In the two-dimensional case of our 
interest this region has an area 7 2 around p + q c . If the 
DOS does not vary appreciably in the Brillouin zone, one 
obtains a result similar to the one found by Engelsberg- 
Schrieffer 



which was derived under the assumption 
of electron-hole symmetry and a q-independent matrix 
element between the electrons and the CM (a Holstein 
phonon): There is a dispersion break at a binding en- 
ergy u> ~ too and a reduction of ojq always leads to a 
shift of the electron dispersion at higher binding energies 
(and only in a limited energy range around the break). 
On the other hand, if p p + qc varies substantially around 
the Fermi surface, a different situation can occur. Let 
us suppose that the isotope replacement acts on 7 [e.g., 
7( ls O) > 7( le O), this is what we find to be the relevant 
effect in our detailed calculations]. Then upon isotopic 
substitution the hole can be scattered in a broader mo- 
mentum region where the average DOS can be substan- 
tially different. In particular, if the local DOS increases 
i?eS( 16 0,w < -loq) < i?e£( ls O,cj < -w ) (cf., e.g., Fig. 
Et) and as a result the electronic dispersion is shifted up- 
wards to lower binding energies Sk( ls O) > Ek( 16 0). This 
is indeed the desired effect to explain the experiments of 
Ref. |3j near the M points. 

In order to substantiate these ideas within more de- 
tailed calculations we consider a system of superconduct- 
ing electrons exposed to an effective action 



S =-9 2 Y] / dn / dT 2 Xq(ri - r 2 )p q (Ti)p_ q (r 2 ) 
Jo Jo 

. (1) 

describing dynamical incommensurate CO fluctuations. 
Using Nambu-Gorkov notation the unperturbed Matsub- 
ara Green's functions G° are 2x2 matrices and the lead- 
ing order one-loop contribution to the self-energy reads 



as 



Xq (ip)T z G° (k - q, «w - ip)r^ (2) 



which in turn allows for the calculation of G = G° + 
G°SG. The spectral function can be extracted from 
ZkR= (l/7r)JmG u (k,w). 

In order to simplify the calculations we consider a 
Kampf-Schrieffer-type model susceptibility which is 
factorized into an u>- and q-dependent part, i.e. 



XqM = W(iw)J(q) 
where, choosing a unit lattice spacing, 
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cos(q x -q%)- cos(q y - q c y ) ' 

(4) 

M is a suitable normalization factor introduced to 
keep the total scattering strength constant while vary- 
ing 7. The (real-frequency)-dependent part W(ui) = 
J dvF^v^vK^-v 2 ) withi^rM ~ T / [{u; - lu ) 2 + 
r 2 ] is a normalized lorentzian distribution function cen- 
tered around ojq with halfwidth T. If only this part were 
present in x(i<jj), one would have a bosonic spectrum 
B(u>) — ta,nh(u>/ (kT))F UOt r(^) which is a "smeared" ver- 
sion of the Holstein phonon considered in Ref. ^3 ■ The 
crucial feature of the susceptibility Q is the substantial 
momentum dependence, which describes the local order 
formation and reflects the proximity to a (missed) CO in- 
stability. This susceptibility contains the charge-charge 
correlations which are enhanced at the four equivalent 
critical wave vectors (±7r/2, 0), (0, ±7r/2). The static 
limit g{v) = 5{v) together with an infinite charge-charge 
correlation length 7—^0 reproduces mean-field results 
for a long-range CO phase ,17]. For the bare electron 
dispersion used in the following 

£fc = — 2i(cos(fe E ) + cos(fc y )) + 4t' cos(k x ) cos(k y ) — u 

the parameters have been fitted so that the resulting 
spectral functions approximately agree in energy with 
the data of Ref. 0. We use t = 0.35eV, t'/t = 0.32, 
u = —0A2eV. The superconducting d-wave gap is 
Afe = 35meV[cos(k x ) — cos(fcj,)]/2. The width of the 
lorentzian in frequency space is T = 50meV 

In Ref. 0] we demonstrated that a static two- 
dimensional eggbox-type charge modulation with q c ori- 
ented along the (1,0), (0,1) directions can account for the 
basic Fermi surface (FS) features in the optimally and un- 
derdoped Bi2212 compounds. In a subsequent paper 
we analyzed ARPES spectra in Bi-2212 showing that the 
main dispersive features 0,0] were reproduced, includ- 
ing the well-known kink in the nodal dispersion which ap- 
pears as a break for momentum scans far from the node 
[19|. The main point to be recollected here is that the 
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M-points define the "hot" regions for incommensurate 
CO scattering. In these regions of k-space the low-energy 
QP's are strongly scattered by the CM's into other states 
with high DOS (still near the M points). On the other 
hand, along T — X the coupling of electronic states to the 
CO mode is significantly reduced ("cold" region) at low 
binding energies and one finds weakly interacting QP's 
which get heavily damped only at high binding energies. 
Here the QP's both those at low and at high binding en- 
ergies are scattered by the CM's into regions with rather 
low DOS. This different character between "hot" and the 
"cold" regions is a key point to understand the behavior 
of the ARPES-IE, with the "surprising" change of sign 
moving from diagonal scans to scans near M and sets the 
stage for our analysis. 

- Only mode-frequency change — First of all we tried 
a simple change of the CM frequency only. Since we 
are not dealing with simple phonons, but rather with 
(nearly critical) CM's, it is difficult to a priori deter- 
mine the isotopic shift (if any) of the CM frequency. We 
take wo( 16 0) = 75meV, wo( 18 0) = 60meV to optimize 
the agreement with experiments. The frequency shift is 
quite large (15meV) with respect to the one expected for 
phonons. The coupling constant is g( le O) = g( ls O) — 
0.l3eV, which leads to self-energies at the nodal point of 
the same order of magnitude than those reported in Q . 
Both along the diagonal direction and along the vertical 
cut corresponding to cut #7 in Ref. Q there is a signif- 
icant isotope shift only in a narrow energy window near 
the kink around — O.leV. Our dispersions do not display 
an isotope shift in the high energy region. Upon chang- 
ing in the opposite way the frequency ivo( 16 0) < ojo( 1s O) 
the agreement with experiments is even worse. There- 
fore, within our model and despite the large shift in loq 
the IE can not be ascribed to a simple shift of the CM 
frequency. 

- Coupling and mode-frequency change — As before 
we choose luo( 16 0) — 75meV and cj ( 18 O) = 60meV. 
In addition we consider now different coupling constants 
for the isotopes, g can be determined from the require- 
ment that the normal contribution to the self-energy 
^u(kF, Ep) at the nodal point gives a reasonable fit to 
the ReT, extracted in Ref. 0. We find g( w O) = 0.13eV 
and g( 1& 0) = O.llAeV. Along the diagonal direction one 
finds qualitatively the same behavior as with luq changes 
only: the IE is still confined to a rather narrow en- 
ergy window. Moreover the theoretical dispersions along 
the vertical cut do not properly describe the findings of 
Ref. 3j: contrary to the experiments, with isotopic re- 
placement the low-energy part of the dispersion for 18 O 
is shifted at higher energies with respect to the disper- 
sion for 16 0. Furthermore, while le O shows the break in 
the dispersion as in the experiment, the lower coupling 
for ls O turns this break into a kink. 

- Only change of 7 — Finally we consider a change 
of 7 alone. We take uj ( w O) = w ( 18 O) = 50meV and 



g( le O) = <?( ls O = 0.16ey, but we choose different spatial 
correlation lengths, i.e. 7 = 0.15 for le O and 7 = 0.30 
for 18 0. Fig. ^shows energy distribution curves (EDC's) 
and the corresponding dispersions for a cut along the 
diagonal (Fig. ^ upper panel) and a vertical cut near 
the M-point similar to cut #7 in Ref. (Fig. flower 
panel, cf. inset to Fig. 01). Along the diagonal direc- 
tion, below the kink energy we find a shift of the 18 O 
dispersion towards larger binding energies. Similarly to 
what is experimentally observed, the IE shift extends far 
below the kink energy. In the low energy regime above 
the kink the shift becomes negligible like in the experi- 
ments. One only finds a minor discrepancy with exper- 
iment since in the crossover region the theoretical 16 O 
dispersion falls below the 18 one. In the vertical cut 
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FIG. I: EDC's and the corresponding EDC-derived disper- 
sions (inset) for the case where only 7 is varied upon isotope 
substitution. Intensities have not been scaled to the same 
amplitudes as done in Ref. 0]. Upper panel: diagonal cut; 
Lower panel: vertical cut at k y — 0.6(7r/a) corresponding to 
cut 7 in Fig. 2 of Ref. (cf. inset to Fig. ^i). 

as shown in Fig. ^ (lower panel) one sees that the ls O 
dispersion is shifted above the 16 O dispersion both at 
low and high binding energies. This behavior is also in 
agreement with the experimental findings of Ref. . We 
notice that, despite the qualitative agreement, the calcu- 
lated shift (~ 15 — 20meV) is substantially smaller than 
the experimental one (~ 30 — 40meV). However, this 
quantitative discrepancy could be solved by assuming a 
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(more realistic) momentum dependence of the QP-CM 
mode coupling g, which instead we take constant for sim- 
plicity reasons. 
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FIG. 2: Real part of the normal contribution to the self-energy 
ReT,u(q,ui). Panels (a,b) show ReT.ii for two momenta along 
the nodal direction and panels (b,c) are for two momenta on 
cut #7 (cf. inset to d). 

Fig. |2 shows i?e£ii(g,w) for the case where 7 is in- 
creased by isotope substitution, while keeping ujq and g 
fixed. Results for two momenta along diagonal cuts are 
reported in the two top panels, where the most prominent 
feature is that at large binding energies the self-energy 
generically decreases upon increasing 7 (i.e. upon de- 
creasing the CO correlation length). Remarkably this 
reduction persists down to very high binding energies. 
While this reduction is expected on general grounds even 
for the phonon case, p|, the frequency extension of this 
effect is rather unusual and arises from the plateau in 
the high-energy part of the self-energy occurring when 
7 is small. On the other hand the self-energy generi- 
cally varies little at low energies, while a moderate in- 
crease, i?e£( 18 0) > ReT,( 16 0), occurs near its maximum 
(around which the kink in the dispersion takes place) due 
to the enhancement of the phase space of the scattered 
QP states induced by the increased 7 ■ This effect is 
even more pronounced in the antinodal region reported 
in the two lower panels. Along cut #7 one therefore ob- 
serves a shift of the 18 O dispersion to smaller binding 
energies. 

- Conclusions — We applied here a simple pertur- 
bative scheme of superconducting quasiparticles coupled 
to charge CM's. From the above systematic analysis we 
find that, if the isotopic substitution changes the CM mo- 
mentum broadening 7, our model accounts both for the 
large energy range in which isotopic shifts are observed 
and for the change in sign of these shifts moving from 
the r — > X cuts to the cuts near the M points. The 
fact that the crucial effect of the isotopic replacement 
is to change the CO coherence length (inversely related 



to 7) also supports the idea that charge fluctuations are 
highly correlated in optimally doped cuprates. This is 
consistent with Ref. where it was found that the 
16 O -^ 18 O replacement shifts the ("missed") quantum- 
critical line to higher dopings (or equivalently confers to 
the system properties of more underdoped materials). 
For optimally or underdoped materials at low temper- 
ature this "pushes" the system more deeply in the nearly 
ordered region away from the "missed" critical line, natu- 
rally reducing the spatial CO correlations. An additional 
broadening of the CM momentum distribution can arise 
from disorder, if the w O — ► 18 (3 substitution is incom- 
plete. Another possible way of arguing can be derived 
from Gutzwiller+RPA calculations of the stripe phase, 
which indicate that upon underdoping the intensity dis- 
tribution of the mode becomes more extended in momen- 
tum space 22] . All these arguments provide support to 
an increase of 7 upon 16 -^ 18 O substitution. 

Since it turns out that isotopic changes of ojq and of 
g play a minor role, one could argue that the CM's not 
necessarily are charge modes (mixed with phonons), but 
could as well be spin modes. Indeed, our perturbative 
scheme previously adopted in Ref. provided an inter- 
pretation of s pec tral features in Bi-2212 quite similar to 
those of Ref. , where the coupling of electrons to an- 
tiferromagnetic (AF) spin CM's was considered. This is 
natural because for both AF and CO scattering the states 
near the M points are the "hot spots", i.e. are most 
strongly affected by the scattering. To (at least prelim- 
inarly) discriminate between charge and spin modes, we 
also checked whether CM's peaked around Qaf = (tt, 7t) 
relevant for AF fluctuations could account for the IE in 
ARPES. While a more detailed analysis of AF modes 
is beyond the scope of this work, we found that CM's 
peaked at Qaf fail in reproducing the major features of 
the experimental isotopic shifts in the dispersion curves. 
Although a more specific analysis for AF spin modes 
should be carried out, our analysis suggests that the crit- 
ical momenta typical of charge modes [stripes and "egg- 
box", with q c ~ (±7r/2,0), (0, ±tt/2)] are an important 
ingredient to account for ARPES experimental findings. 
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